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The Genesis of an Impulsive Coronal Mass Ejection observed at Ultra-High 

Cadence by AIA on SDO 
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ABSTRACT 

The study of fast, eruptive events in the low solar corona is one of the science ob- 
jectives of the Atmospheric Imaging Assembly (AIA) imagers on the recently launched 
Solar Dynamics Observatory (SDO), which take full disk images in ten wavelengths 
with arcsecond resolution and 12 sec cadence. We study with AIA the formation of 
an impulsive coronal mass ejection (CME) which occurred on June 13, 2010 and was 
associated with an M1.0 class flare. Specifically, we analyze the formation of the CME 
EUV bubble and its initial dynamics and thermal evolution in the low corona using 
AIA images in three wavelengths (17lA, 193A, and 21lA). We derive the first ultra- 
high cadence measurements of the temporal evolution of the CME bubble aspect ratio 
(=bubble-height/bubble-radius). Our main result is that the CME formation under- 
goes three phases: it starts with a slow self-similar expansion followed by a fast but 
short-lived (~ 70 sec) period of strong lateral over-expansion which essentially creates 
the CME. Then the CME undergoes another phase of self-similar expansion until it 
exits the AIA field of view. During the studied interval, the CME height-time profile 
shows a strong, short-lived, acceleration followed by deceleration. The lateral overex- 
pansion phase coincides with the deceleration phase. The impulsive flare heating and 
CME acceleration are closely coupled. However, the lateral overexpansion of the CME 
occurs during the declining phase and is therefore not linked to flare reconnection. In 
addition, the multi-thermal analysis of the bubble does not show significant evidence of 
temperature change. 

Subject headings: Sun: coronal mass ejections (CMEs) 

1. Introduction 

Synoptic observations of Coronal Mass Ejections (CMEs) have recorded thousands of events 
over the last 15 years leading to a decent understanding of their basic physical parameters (e.g., 
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Schwenn et al.l (|2006l )). However, their initiation and driving mechanisms remain unclear (jForbes et al 
2000 ) because their formation and initial stages of evolution occur in the low corona where, until 
recently, observational coverage was sparse and cadence was low. 

Observationally, we can decompose the CME evolution into three stages: initiation, acceler- 
ation, and propagation (jZhang &; Derd l2006l ). Height-time measurements of CME fronts, using 
EUV imagers and white light coronagraphs, have shown a close correspondence between the initial 
CME acceleration phase a nd the impulsive phase of the soft X-ray (SXR) flare (jZhang et al.ll2004l ; 
Temmer et al.l 120081 . l2010l ) and have established the CME acceleration profile as a proxy to the 
energy release. Detailed acceleration profiles, especially for impulsive CMEs, provide important 
constraints for all CME initiation models. 

The understanding of CME formation is lacking behind the understanding of CME acceleration 
due to cadence and field of view restrictions. Most CME models agree that the final ejected structure 
is a ma gnetic fluxrope which may corres pond to the cavity observed in 3-part CMEs in the outer 
corona (jSubramanian fc Vourlidas! 120071 ) . But we do not yet know if the fluxrope preexists, if it 
forms on-the-fly or if the final structure is a combination of the above processes. The answer 
holds serious implications for the relative importance of ideal versus non-ideal mechanisms for 
the eruption process (e.g.. iDemoulin Sz Aulanierl 120101 . and references therein). Impulsive CMEs 
are of particular interest because they tend to be associated both with strong flares and particle 
acceleration. However, the formation of impulsive CMEs is a particular challenge for existing 
instrumentation. They form low in the corona amidst significant pre-existing fine scale structure, 
they interact violently with their surroundings generating disturbances (such as waves and loop 
deflections that can obscure the evolution of the actual CME), and most of their energy release 
occurs within very short time intervals, at or below the typical cadence of most EUV or white light 
imagers. 



Patsourakos et al.l (J201Cj) (hereafter, PVK10) analyzed relatively high cadence (75 - 150 sec) 



EUV STEREO observations of an impulsive CME. The CME had the appearance of a bubble and 
it was shown that the bubble was an early instance of the CME fluxrope seen higher up in the 
coronagraph field of view. The wide separation between the STEREO observations (~ 48°) allowed 
PVK10 to perform detailed 3D modeling of the erupting bubble. The modeling showed that the 
bubble underwent a phase of lateral overexpansion (i.e., the bubble radius increasing faster than 
its height) very early in its evolution which coincided with the impulsive phase of the flare. The 
overexpansion transformed the bubble from a small, arcminute-scale feature into a large CME- 
scale fluxrope structure within 5 minutes. If this result holds in other events, it may explain the 
discrepancy between the sizes of CMEs and post-eruptive loops arcades and will clarify the true 
relationship between CMEs and their source regions. To accurately identify, delineate and measure 
the size of the erupting bubble, and the subsequent fluxrope, we need high cadence, high sensitivity 
off-limb EUV images which only the EUV imagers onboard the Solar Dynamics Observatory (SDO) 
can provide. In this Letter, we present the first multi-wavelength analysis of the early stages of a 
very impulsive CME with ultra-high cadence (12 sec) and compare the 3D evolution of the erupting 
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bubble to the impulsive phase of the flare. Using three of the ten AIA channels sensitive to different 
temperatures we search for evidence of plasma heating. The observations and data analysis are 
presented in Section 2, and the results are discussed in Section 3. 



2. Observations and Data Analysis 



We use EUV imaging observations from AIA (AIA; iTitle et al.l 12006!) on board the recently 
launched SDO and the Extreme Ultraviolet Imager (EUVI; I Wuelser et al.ll2004l ) onboard the Solar 
Terrestrial Relations Observatory (STEREO) Ahead spacecraft (STA). The eruption occurred on 
June 13, 2010 in active region NOAA 11079 which was located slightly behind the solar limb 
(S24W91). The eruption was accompanied by a rather short duration M1.0 flare starting at 
05:33 UT, peaking at 05:39 UT and returning to pre-event levels at around 06:43 UT. Associ- 
ated with this eruption was a small CME seen by the STEREO coronagraphs. STA was situated 
74° west from Earth, so AR11079 was 17° east of STA's central meridian. The region was therefore 
observed in quadrature by the AIA and EUVI telescopes. 

The eruption was observed as a limb event in all AIA channels with 12 sec cadence but here 
we use only the coronal channels centered at 17lA, 193A and 21lA (Level 1.5 images) and the 
EUVI-A 195 A channel (henceforth, we refer to a given channel simply by its central wavelength). 
The EUVI-A 195 cadence was 2.5 min. 

The active region exhibited considerable activity throughout the day with loops rising and 
brightening but no major ejection was observed before our event. We focus on the eruption, and 
concentrate on the time interval from 05:30 to 05:43 UT. We provide two movies online in AIA 193 
and EUVI-A 195 (moviel.mpg, movie2.mpg). Snapshots fro m these movies are s hown in Figure [TJ 



The image contrast is enhanced with a wavelet technique (jStenborg et al.ll2008l ). Only the large 
scale brightness component has been removed. 

The EUVI-A and AIA observations show the slow rise of loops over the active region for several 
minutes before the CME. Starting at ~ 05:33 UT, brightenings appear very close to the surface 
at the inner core of the region. At the same time, a filament is clearly seen rising and heating 
(Figure [H middle column). The filament seems to rotate as it rises which may indicate that it 
is kinking or writhing. The rise of the filament is associated with the rise of surrounding loops 
which eventually form an EUV cavity resembling a bubble around 05:36 UT. Taking advantage of 
the two viewpoints afforded by the EUVI-A and AIA observations we use triangulations to locate 
the erupting filament and the bubble in the two telescopes (see features marked in the righttmost 
column of Figure [T]) . 
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2.1. Formation of the bubble 

Thanks to the high cadence of the AIA images, we can study the formation and evolution 
of the EUV bubble in great detail. We provide online a base-ratio movie in 171 (movie3.mpg), 
where we divide each wavelet-enhanced image by a pre-event image at 05:30:23 UT. Representative 
snapshots from the movie are in Figure [2j The motion of the rising filament appears as twisting 
motion in AIA an d straightening m otion in EUVI-A. This pattern seems consistent with writhe 



rather than kink (jTorok et al.ll2010i ). As seen in AIA, the filament rises, twists and then starts 
falling towards the surface (after ~ 05:33:11 UT). At that time, surrounding loops with seemingly 
different orientations begin to slowly grow (e.g., alternating black- white patterns in the lower two 
panels in Figure [2|). Simultaneously, a ribbon-like intensity front propagates away from the filament 
along a south-east direction (upper two frames, Figure [2|). The extent of this front roughly matches 
the lateral extension of the expanding loops. In the meantime, the expanding loops seem to pile 
up at the surface of the forming bubble (upper panel of Figure [2]). The bubble becomes clear at 
~ 05:36 UT and continues to grow. After 2 minutes the bubble begins a strong lateral expansion 
which we analyze later. This time also signals the launch of an EUV wave around the bubble 
(moviel.mpg) but the wave analysis will be reported elsewhere. The outer rim of the bubble 
becomes progressively thinner and dimmer after ~05:38:23 UT. After ~05:42 UT, it increasingly 
becomes difficult to trace its upper part. The bubble exits the AIA field of view at ~05:45 UT. 



2.2. Evolution of the bubble aspect ratio 

Following PVK10, we follow the evolution of the bubble aspect ratio to quantify the early 
3D evolution of the CME. Unfortunately, the very impulsive nature of the event does not permit 
detailed two- viewpoint fitting of the bubble, as was done by PVK10. There is only one frame with 
the bubble fully formed in both EUVI-A 195 and AIA 193 (third column, Figured}. 

We are constrained to use the AIA data only. Our approach is to fit circles to the bubble 
rim using the 171 base-ratio images for the interval 05:35:23-05:42:23 UT. There are 36 images 
in all. We do not use other times because it was difficult to identify a significant fraction of the 
bubble rim with reasonable confidence due to either its slow evolution or faint signal. For each 171 
base-ratio image, we manually select 10 points along the bubble rim and fit a circle to these points. 
We assume an error of 5 AIA pixels (~ 3 arcsec) in the determination of each point. This rather 
conservative estimate is the median of the half- width of the bubble rim for the 36 images. Figure [3] 
shows examples of the circle fits and demonstrates that a circle represents a good approximation 
to the upper section of the bubble. To further check the consistency of the circular assumption, 
we perform a circular fit to the EUVI-A bubble at 05:38:00 UT which yields a radius within 20% 
of the AIA 193 fit at 05:38:06 UT (Figure [I]). This agreement indicates that the bubble can be 
reasonably approximated by a 3D sphere. 

The circular fits supply the radius and center height, relative to the solar surface, of the bubble. 
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Hence the ratio between these two quantities provides the aspect ratio of the erupting bubble 
(Figure Ufa)). The error bars are determined by standard error propagation of the uncertainties 
of the height and radius measurements. If we assume that the bubble represents a rising fluxrope 
seen edge-on, we can interpret the derived aspect ratio as the ratio of the major to the minor axis 
of the erupting fluxrope. Q| 

The aspect ratio has a constant value of ~ 1.7 between ~ 05:36:00 UT and ~ 05:38:20 UT. 
The ratio decreases rapidly, within approximatively 70 sec, to a value of ~ 1 and remains relatively 
constant thereafter. The "jump" in the aspect ratio before ~ 05:36 UT is possibly due to the 
uncertainty in the bubble determination at that time. 



2.3. Bubble acceleration profile and flare dynamics 



Following common practice, we measure the height-time (HT) evolution of the EUV bubble 
in 171 by tracking a point in the bubble front between 05:34:00-05:43:11 UT. The resulting 47 
measurements are shown in Figure \Mb). As done freq uently, the HT data are first smoothed to 
reduce small-scale fluctuations (e.g.. TVrsnak et aL 2007 ). We use a smoothing cubic spline scheme 
which minimizes a function consisting of the sum of a x 2 fit of the data with a cubic spline plus 
a penalty function proportional t o the second der ivative of the cubic spline multiplied by a user- 
supplied factor < spar < l,e.g. IWeisberd (]2005l ). Smoothing factors between 0.7 — 0.8 provide 
the best compromise between noisy (zero or small spar) and extremely smoothed (large spar) 
acceleration profiles. We adopt the median of the above interval (spar = 0.75) for the smoothing 
of the HT data (Figure 0(b)) which are then differentiated with respect to time to obtain velocity 
and acceleration profiles (panels (c) and (d) of Figure 0}. The error bars are derived from the 
maximum residuals between the curves calculated with the median spar and those for spar 0.7 
and 0.8. The errors in the speed and acceleration profiles resulting from the use of 5-pixel error 
estimate in the HT measurements were very small and were thus neglected. The results show 
that the bubble accelerates rapidly (within 4 min) to a speed of ~ 400 kms -1 and then starts 
to decelerate (Figure 0b). The acceleration profile is characterized by a short-lived (FWHM ~2 
min) acceleration pulse followed by deceleration. The latter is e xpect ed to occur even in the inner 
corona, for impulsively accelerated CMEs (e.g., iTemmer et al.l (|2010l )). This may be due to solar 
wind drag as the CME speed exceeds the local solar wind speed. The deceleration is verified by 
the LASCO C2 measurements which show a speed of ~ 200 kms' 1 . 

The acceleration profile corresponds closely to the time derivative of the GOES soft X-ray 
(SXR) light curve, a proxy for hard X-rays (HXR) and hence for the flare r econn ectiorH, (Fig- 
ure H{d)), as generally expected (e.g., IZhang et al.l 12004 : ITemmer et al.l 120081 . l2010i ). The slower 



lr The spherical approximation applies only to the upper section of the flux-rope and not to its legs. 
2 RHESS1 was observing the Crab Nebula during our event. 
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SXR rise profile arises from the occultation of a significant part of the flaring region from Earth. 
The SXR loops contribute to the GOES lightcurve only after they rise above the solar limb. The 
SXR levels are high when the bubble forms (~05:36 UT). Both the SXR derivative (and hence re- 
connection) and the bubble acceleration are close to zero when the strong lateral expansion starts. 
It is an indication that the expansion may be driven by mostly an ideal process. We will return to 
this important point in the Discussion. 



2.4. Multi-temperature behavior of the bubble 

We investigate the thermal behavior of the expanding bubble using the almost simultaneous 
AIA observations with filters sensitive to coronal plasmas of different temperatures. We use the 
observations in 171, 193, and 211 channels with peak responses at 0.8, 1.25, and 1.6 MK, respec- 
tively. 

We posted an online movie of the 171/193 intensity ratio during the event (movie4.mpg). 
The observations between the two channels are separated by only 5 sec which corresponds to a 
displacement of around 4 pixels at the maximum speed of ~ 300 km s" 1 . This separation is smaller 
than the 10 pixel full width of the erupting rim and so it has minimal impact in the calculated 
intensity ratios. Snapshots from the intensity ratio movie are shown in Figure [5l Brighter (darker) 
areas could signify relatively larger (lower) amounts of warmer (hotter) plasmas Q Such areas seem 
to remain unaltered which implies that the temperature does not vary substantially during the 
event. In Figure [5j the bright/dark segments of the bubble remain as such during its evolution. 

To further quantify the thermal evolution of the bubble, we construct lightcurves by extract- 
ing the intensity along the propagation path between 05:30 UT and 05:42 UT in each channeQ. 
The intensities at 5:30 UT are subtracted from each exposure. The peak intensity for the base- 
subtracted lightcurves is used as a proxy for the bubble intensity. All three lightcurves peak nearly 
simultaneously and at about the same time as the strong lateral expansion of the bubble (Figure H]). 
Also, the lightcurves have similar lifetimes to the acceleration pulse (FWHM ~ 2 min). The 171 
lightcurve, however, increases at a slower rate than either the 193 or 211 curves. This may be 
because the loops rising from the AR core were emitting mostly in 171. The hotter emissions seem 
to lie around these 171 loops and the bubble in those channels seems to form as a result of the 
initial expansion of the 171 structures. The discussion above suggests that the bubble temperature 
did not vary substantialy during its expansion. 



3 Note that filter-ratios are not strictly monotonic functions of temperature over large temperature intervals. 

4 The intensities in the hotter channels, i.e. 335, 94, 131 were affected by diffraction patterns from the flare emission 
and were unsuitable for quantitative analysis. 
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Discussion and Conclusions 



Arguably, our most interesting finding is the strong lateral overexpansion of the bubble that in 
essence 'inflates' an initially small-scale feature into a CME-scale structure. This period of inflation 
starts after the main phase of flare heating and CME acceleration and argues strongly for ideal MHD 
rather than a reconnection-related process. The observed overexpansion is a considerable challenge 
for any CME model. The overexpansion may be triggered when the expanding CME 'exits' above 
the AR loop core and encounters weaker magnetic fiel ds. To achieve pr essure balance with its 
surroundings, the CME bubble must expand laterally (jMoore et al.l 120071 ). Another explanation 
may be flux conservation around a rising fluxrope of decreasing current (PVK10). However, the 
overexpansion in PVK10 was relatively well-synchronized with the impulsive flare phase and hence 
a reconnection-based driver was possible. The differences between these impulsive CME events 
indicate that the characteristics of lateral overexpansion in the early phase of CMEs vary from 
event to event making their study a potentially sensitive diagnostic of CME evolution and coronal 
properties. For instance, a deeper study of this event may supply estimates on the magnetic field 
strength of the regions in the vicinity of the eruption site. A detailed MHD modeling of the PKV10 
event to determine the relative role of ideal vs non-ideal processes is underway. 

Our multi-temperature observations show no appreciable temperature changes in the bubble 
during the event. If significant plasma heating (or cooling) was taking place, it would ha ve led 
to temporally displaced lightcurves in the different channels (e.g., iRobbrecht fc Wanal2010l ). This 
does not necessarily preclude heating (and then cooling) of the bubble plasma during its forma- 
tion. According to standard flare-CME theory, field lines reconnecting i n the vertical cu rrent sheet 
underneath the erupting flux become part of the CME bubble (e.g., iLin et al.1 12004 ). If these 
lines fill with hot plasma (i.e., by evaporation) at a slower rate than the eruption, then no tem- 
perature changes may be detected. Indeed, for the analyzed AIA chan nels (« 0.5-3 MK), the 
evaporation speeds are around 100 kms -1 (e.g.. iMilligan Sz Dennis! 120091 ) — slower than the CME 
speed(~ 200 — 400 kms -1 ). More analysis is needed to test these ideas further. 

While reconnection may not be important in the lateral overexpansion of the bubble, it could 
be a significant factor for its formation. The bubble forms from a set of pre-existing loops at 
varying orientations during the main flare phase. This is consistent with the trans formation of 
a loo p arcade to a fluxrope structure and is predicted by several CME models (e.g., iLynch et al 



20081) . However, the possible eruption trigger lies in the rise and possible instability of a small 
filament. Such structure could be related to a pre-existing fluxrope. At any rate, our observations 
of the initially very small bubble radius (~ 0.03-Rq) set an upper limit for the size of any pre-existing 
flux rope. We suggest that the observed lateral overexpansion may well the process through which 
eruptions starting small in the corona become large-scale CMEs further out. Analysis of more 
events is required to establish or refute this suggestion. 

To summarize, we have presented EUV observations of the genesis of an impulsive CME. By 
taking full advantage of the unique SDO AIA capabilities — its unprecedented high cadence (12 
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sec) and multi-temperature coverage — we were able to resolve the various stages of the event in 
great detail. Our main findings are: 

• A set of slowly rising loops, possibly triggered by a rising and maybe kinking or writhing 
filament, is transformed into an EUV bubble within 2 minutes. 

• The EUV bubble forms when both flare heating and CME acceleration are at their maximum 
levels. 

• The bubble experiences a 70 sec period of strong lateral overexpansion followed by self-similar 
evolution. 

• The lateral over-expansion starts when flare reconnection and CME acceleration are well 
through their peaks. 

• The bubble rim emission shows no significant evidence of temperature change. 

The AIA data used here are courtesy of SDO (NASA) and the AIA consortium. We thank the 
AIA team for the easy access to calibrated data and the referee for useful comments. The SECCHI 
data are courtesy of STEREO and the SECCHI consortium. 
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Fig. 1. — EUV Observations of the erupting active region in AIA 193 (upper panel) and the EUVI- 
A 195 channel (lower panel). An untwisting filament apparently associated with the initiation of 
the eruption is marked by arrows. The dashed lines outline the expanding bubble. EUVI-A is 
situated 74° West from Earth. 
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Fig. 2. — Representative AIA base-ratio 171 frames before and towards the EUV bubble formation. 
The ratio increases with from black to white within the range 0.3-1.4. 
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Fig. 3. — Representative circular fits (green lines) of the bubble in AIA 171 base-ratio images. The 
blue box marks the bubble center. The image at 05:30:23 UT is used as the base. 
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Fig. 4. — (a) Temporal evolution of the aspect ratio of the bubble, (b) HT measurements of the 
bubble front (crosses). The points are connected with spline smoothing curve (solid line), (c) Speed 
versus time curve prouced by differentiating the smoothed HT measurements, (d) Acceleration 
versus time curve (crosses) over-plotted on the time-derivative of the GOES SXR curve (solid line). 




Fig. 5. — (a)-(c) AIA 171/193 intensity ratio snapshots, (d) Light-curves across the bubble rim in 
171, 193, and 211 with pre-event intensity levels subtracted. The times on the images correspond 
to 171 and the 193 images trail by 5 sec. 



